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Abstract 

Barrientos Bendezu and Kniehl |J recently suggested that W^H^ associated 
production may be a useful channel in the search for the elusive heavy charged 
Higgs bosons of the 2 Higgs Doublet Model at the Large Hadron Collider. We 
investigate the phenomenology of this mechanism in the Minimal Supersym- 
metric Standard Model, with special attention paid to the most likely heavy 
Higgs decay, H T — > tb — > bbW^, and to the irreducible background from 
top pair production. We find that the semi-leptonic signature l bbW + W~ — > 
bb jj t + missing momentum' is dominated by top-antitop events, which over- 
whelm the charged Higgs signal over the heavy mass range that can be probed 
at the CERN collider. 
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1. Introduction 



In Ref. m, Barrientos Bendezii and Kniehl pointed out that the processes 

bb -> W ± H T (1) 

and 

^ - W^if* (2) 

can be used in the search for the heavy (M#± > m t + mj) charged Higgs bosons H ± of 
the 2 Higgs Doublet Model (2HDM) at the Large Hadron Collider (LHC). 

This result is particularly welcome, since it has been remarked in several occasions 
(see Ref. for an overview) that at present it is not at all certain that such particles 
can be detected at the LHC, even in the Minimal Supersymmetric Standard Model 
(MSSM), if the typical energy scale of Supersymmetry (SUSY) is much greater than 
the charged Higgs mass. 

Previous studies have shown that, if M H ± > m t + m&, the chances of H ± detection 
at the LHC are reliant only on two production mechanisms: the subprocesses bg — > tH ± 
U and bq -> bq'H* [g and provided that M H ± Z 300 - 400 GeV f|. These channels 
generally have poor signal-to-background ratios, as the event signatures always involve 
a large number of jets, which is the typical noise of a hadron-hadron machine. 

When H ± 's are heavy, for Msusy 3> M H ±, they decay almost exclusively to bt{bt) 
||. In addition, hadronic signatures of W ± bosons produced in top decays are normally 
selected in order to allow for the H ± mass reconstruction. Therefore, at first sight, it 
appears that the signals ([I])-(|2]) advocated in Ref. [|IJ as a useful source of Higgs events 
may be swamped by the irreducible background from top pair production, 

qq — > ti and gg — > tt, (3) 

with subsequent decay through the intermediate statef] bbW + W~ . In order to under- 
stand whether this is the case, we made a detailed signal-to-background analysis and 
found that top-antitop production and decay indeed overwhelms the new Higgs signal. 

The plan of this paper is as follows. In the next Section we describe the details 
of the calculation and list the parameter values adopted. Section 3 is devoted to the 
discussion of results. We present our conclusions in Section 4. 

2. Calculation 

We generated the signal cross sections by using the formulae of Ref. QTJ . However, for 
phenomenological analyses we need to supplement those expressions in several ways. 

Firstly, their matrix elements (MEs) do not include the boson decay and thus 
carry no information on the angular distributions of the fermion pair it produces. For 
the this does not matter, since scalar particles decay isotropically. However, even 

1 Alternatively, in a narrow window in M H ± and only at low tan/3, the charged Higgs bosons can 
decay to W^h pairs, where h represents the lightest Higgs boson. Although we will not treat this case 
here, we note that even in this channel the final state is identical to that of top-antitop, as h — » bb is 
dominant over most of the SUSY parameter space. 



in this case one has to provide the correct kinematics for the decay sequence H — > 
tb — > bbW^ 1 — > bbjj (where j represents a light quark jet produced in the W ± decay), 
which we have done by computing the exact ME constructed by means of the HELAS 
subroutines 0. As for the W ± decay in the production channels (0) and (§), we have 
re-evaluated their MEs with the additional insertion of the W boson decay currents. 
For the bb fusion case, the actual expression is the same as for the bq — ► bq'H ± process 
calculated in @] and recently modified in |§, but with some leg crossings. For the gg 
fusion case, the result is simply the replacement 

\(s,M^ ± ,M 2 H± ) -> VM^ ± |^±| 2 [2(pi -Ph±)(P2 -Ph±) ~ (pi ■P2)M 2 H± ] (4) 

in equation (8) of 0, where g 2 = 47ra em /sin0^, \G W ±\ 2 = {{Pw± - M^ ± ) 2 + 
{Yw ± ^w ± ) 2 Y 1 1 with Ph±, Pw ± , Pi an d Pi the four-momentum of the H ± , W ± , first 
and second lepton, say £ and U£, from the gauge boson decay, respectively. 

Secondly, their MEs for the gg fusion processes do not involve squark loops, this 
preventing one from studying possible effects of the SUSY partners of ordinary quarks, 
when their mass is below the TeV scale. For example, these corrections are expected 
to be sizable in the MSSM, which is adopted here as the default SUSY framework. 
In this respect, we have modified the gg triangle formula of Ref. for the case of 
intermediate neutral Higgs production (top graph in Fig. 2 there), by inserting the 
well known H [| |l(J squark loop terms p. 

Conversely, we have not included here the contribution of the box diagrams (bottom 
two graphs of Fig. 2 in Ref. [l|]), for which the authors of that paper found no compact 
expression. According to their curves in Fig. 5 [|TJ , this should result in an overestimate 
of the total cross section of subprocess (0), as the triangle and box diagrams interfere 
destructively. The overall effect is however negligible at large tan/3, the regime where 
the W ± H Zf cross section is larger (whereas, for tan/5 = 1.5, it can at times be more 
than a factor of two: see Fig. 3 of Ref. |l|]). In addition, this is particularly true for 
Mh± 3> m t , the mass interval with which this paper is concerned. Indeed, in most of 
our plots we will concentrate on that portion of the MSSM parameter space, for which 
we expect our results to be reliable. 

Finally, all our calculations for the signal were tested against the original cross 
sections of [JJ, and also using MadGraph |TTJ for the case bb — > W^H^. 



For the background we have assumed that the QCD contribution is reducible by 
cuts on the reconstructed top and W ± masses. Therefore, we studied the top pair 
background only, which is anyhow the dominant component of the final state bbUD£i?£ 
where U and D refer to up- and down- type massless quarks and I = e or \x. We have 
considered both qq —>■ tt —>■ bbVDivi and gg — > tt — > bbUDiPg, which we generated 
at leading order using the HELAS library 0. The outputs of the corresponding code 
agree with the results given in Ref. |12[ when the W s are on the mass shell. In fact, 



notice that finite widths effects of top quarks, gauge and charged Higgs bosons have 
been taken into account here. 

Concerning the values of the various MSSM parameters entering the computation 
of the signal processes ([TJ)-([2]) (and, marginally, the MSSM top width), we proceeded 
as follows. First, we produced the masses and the widths of the Higgs bosons by means 
of the two loop relations of Ref. (^3] (see also @). To simplify the discussion, we have 



assumed a universal soft Supersymmetry-breaking mass m~ = m"~ = m~, and negligible 



mixing in the stop and sbottom mass matrices, A t = Ab = fi = 0. Second, the MSSM 



Higgs widths were generated using the program HDECAY |L5], which in turn uses mass 
relations at the same perturbative level. Squark masses entering the loops in the gg 
induced signal processes have been kept as independent parameters and their values 
varied between 300 GeV and 1 TeV, the minimum figure being chosen in such a way 
that the superpartners do not enter the H ± decay chain, as for the upper mass of the 
latter we have taken the value of 600 GeV. Further notice that squark masses have been 
considered degenerate, for illustrative purposes, so that only sbottom and stop loops 
in practice contribute. 

In the numerical calculations presented in the next Section we have adopted the 
following values for the electromagnetic coupling constant and the weak mixing angle: 
a em = 1/128 and sin 2 ^ = 0.2320. The strong coupling constant a s , which appears in 
higher orders in the computation of the charged Higgs decay widths and enters in some 
of the production mechanisms, has been evaluated at one or two loops, depending on 
the Parton Distributions Functions (PDFs) used, with (for the number of active 
flavours Nf = 4) input according to the values adopted in the fits of the latter. For the 
leading-order (LO) package CTEQ4L [fTEjl , which constitutes our default set (as in [|TJ), 
we have taken 236 MeV. The factorisation/renormalisation scale Q entering in both a s 
and the PDFs was set to Vs, the centre-of-mass (CM) energy at the partonic level, in 
all processes generated. 

In order to have an estimate of the dependence of the bottom quark and gluon 
structure function we tested our signal rates against the next-to-leading (NLO) sets 
MRS(R1,R2,R3,R4) |T7|, i.e., the Martin-Roberts-Stirling packages of the same 'gen- 



eration' as the CTEQ ones considered here, plus the newly presented sets MRST 
which embody new data and an improved description of the gluons at small x, along 
with a dedicated treatment of the heavy quark structure functions. Typical differences 
were found to be within 15-20%, at cross section level, with the shape of the relevant 
differential distributions being little affected by the treatment of the partons inside the 
proton. 

For the gauge boson masses and widths we have taken Mz = 91.1888 GeV, Y z = 2.5 
GeV, M w ± = 80.23 GeV and T w ± = 2.08 GeV. For the top mass we have used 
m t = 175 GeV with the corresponding width evaluated at tree- level in the MSSM 
(yielding T t = 1.55 GeV if M H ± > m t — nib, the Standard Model value). Bottom 
quarks have been considered massless when treated as partons inside the proton, while 
a finite value of 4.25 GeV (pole mass) has been retained in the final states. Note 
that in calculating the ME for the decay process if ± — > tb the Yukawa mass of the b 
quark has been run up to the charged Higgs mass scale, in accordance to the way the 
corresponding width has been computed. Finally, for simplicity, we set the Cabibbo- 
Kobayashi-Maskawa matrix element V^km to one. 



3. Results 

As it is impractical to cover all possible regions of the MSSM parameter space 
(Ma, tan/5), we concentrate here on the two representative (and extreme) values 
tan/3 = 1.5 and 30, and on masses of the charged Higgs boson in the range 160 GeV 



~ Mh± ~ 600 GeV. The large bibliography existing on the MSSM Higgs decay phe- 
nomenology should allow one to easily extrapolate our results to other values of tan f3 

We consider the semi-leptonic event modes for processes (ID)-©: 

W ± H Zf — > W^tb — > bbW + W~~ — » bb jj £ + missing energy/momentum. (5) 
In fact, we base our signal selection procedure on the following general requirements. 

1. High purity double b quark tagging. This may be expected to yield an efficiency 
of at least 50% per fiducial b jet |19|, pOf] . This is essential considering the large 
rate of W ± + jet events with light quarks and gluons. All our results will assume 
100% bottom quark tagging efficiency e&, and thus will eventually need to be 
multiplied by the actual el once we will have the LHC detectors running. 

2. Lepton isolation at high transverse momentum. Selecting semi-leptonic events 
should enable one to use the high and isolated lepton (electron and/or muon) 
originating from the W T produced in association with the H T as a clean trigger^ 

In addition, the light quark jets j coming from the secondary W T , from 
H T — > tb — > bbW T , would allow for the reconstruction of the Higgs mass peak. 
However, one should recall that the two gauge bosons could decay the other 
way round, this in principle spoiling the efficiency of the signal selection. In 
practice, one can remove the contribution from semi-leptonic Higgs decays by 
simply imposing cuts on the reconstructed top mass. 

3. W ± and t mass reconstruction in two and three jet combinations, respectively, to 
eliminate QCD multi-jet production. 

Our results are shown throughout Figs. [1|-|8|. When differential spectra are plotted, 
the three representative parameter space points of Ma = 200, 400 and 600 GeV at 
tan /3 = 30 have been chosen, corresponding to M#± = 214, 407 and 605 GeV, respec- 
tively. 

As a preliminary exercise, in order to understand the kinematics of the signal and 
background better, and possibly to pin down systematic differences which can be used 
in the selection of candidate Higgs events, we compare their total and differential rates 
in the channel fl5|) without the usual detector cuts on transverse momenta and pseudo- 
rapidity. In fact, both processes have finite production rates at lowest order over all of 
the phase space. 

The main frame of Fig. |T] shows the signal (with no squark loop contributions) and 
background cross sections before any detector or selection cuts, evaluated at the LHC 
energy (14 TeV) and plotted against the H ± mass for the two values of tan/?. The 
background remains constant as a function of the charged Higgs mass, except when the 
decay t — ► bH + is kinematically allowed. It is found that even when the signal cross 
section is the highest the signal/background ratio is less than one in a thousand. Notice 
that at decay level, see eq. @, the signal rates suffer from a further small depletion 
(in addition to that due to finite width effects) as compared to the on-shell ones of 

2 For the time being, we neglect W — > tv t decays, which should also be identified easily thanks to 
their 'one-prong' signatures, as remarked in H. 



Ref. PJ, if the latter are multiplied by the relevant branching ratios (BRs). This is due 
to the fact that the Goldstone part of the longitudinal component of the W ± boson 
(i.e., the wave- function term proportional to p^ v ±/M w ±) produced in association with 
the if T scalar in hh fusion (included in that reference) does not survive the decay into 
'massless' fermions. The effect depends on tan/3, owning the structure of the MSSM 
couplings in the ^W ± H Zf vertices, where $ = H and h, and on the relative strength 
of the two contributing diagrams (see Figs. 1-2 of Ref. [0]). 

In the central insert of Fig. |l] we study possible virtual effects of SUSY, manifesting 
itself in the triangle diagrams of gg fusion. In fact, we plot the ratio of the signal 
cross sections obtained by adding the rates of both suprocesses (H) and (§), the latter 
including squark loops, against those calculated when such contributions are neglected. 



As already remarked in the literature fllQ"|, sizable SUSY effects in the gg subprocess 



are expected only for squark masses below 500 GeV or so and particularly at small 
tan p. Thus, in our plot we present the ratios for rrig = 300 and 500 GeV at tan (3 = 1.5 
only. Given the remarks made in the Introduction, concerning our remotion of the box 
diagrams and, consequently, of the cancellations against the triangle ones (which we 
would further think to be active for virtual squark contributions as well), our numbers 
should in this circumstance be intended as a sort of upper limit that one might expect 
from such SUSY effects. From this prospect, it is then clear that production rates 
can increase by no more than 20% or so, and provided squark masses are rather low, 
a correction indeed comparable to the combined uncertainties related to the b quark 
and gluon PDFs in (P and @, respectively, and to the scale dependence of the gg 
production rates jlTJ. For this reasons, and to simplify the discussion as well, hereafter, 
we will neglect altogether the squark loop contributions in our signal rates. 
The differential spectra are displayed in Figs. |2|-|5]. In particular, we plot[| 

• the lepton, light quark, b jet and missing transverse momentum; 

• the lepton, light quark and b jet pseudorapidity; 

• the lepton/light quark jet and lepton/6 jet separation, defined by the variable 



AR = y (A77) 2 + (A0) 2 in terms of pseudorapidity 77 and azimuth 0. 
Furthermore, we present the invariant mass spectra of the following systems: 

• bb, as obtained by pairing the two jets with displaced vertices; 

• three jets, with only one b jet involved; 

• four jets, involving all four jets in the final state. 

The combinatorics in the three, four and lepton/light quark jet systems is accounted 
for by simply plotting all possible momentum combinations each with the same event 
weight. In other terms, the signal spectra contain both hadronic and leptonic W ± /H Zf 
decay modes (but not their interference, that we expect negligible), which have in fact 
different kinematics. The two component will eventually be separated by imposing 

3 We make no distinction between b and b jets, tacitly assuming neither jet charge determination 
nor lepton tag. 



cuts around the reconstructed top quark mass, as mentioned before. In contrast, the 
background spectra, obviously, do not depend on whether the top or the antitop decays 
leptonically. 

The usual distributions in p? and rj indicate the effects of detector acceptance cuts on 
the signal and background samples. Neither of these affects the event rates significantly. 
The distribution in AR indicates that the requirement of lepton/jet separation will not 
harm the event rates either. Furthermore, the signal and background distributions 
are very similar and it is clear that none of these variables can profitably be used to 
optimise the selection procedure. Incidentally, we mention that we also had a look at 
the transverse momentum and pseudorapidity of the three hadronic systems introduced 
above, without finding any significant difference between Higgs and top events. 

Presumably, the invariant mass distributions, see Fig. [5] and ^ for the purely hadronic 
and semi-leptonic systems, respectively, will give us the greatest chance of removing 
the background. By imposing cuts on the two light quark jet and two light quark plus 
bottom jet invariant masses around the and top quark resonances, respectively, we 
can remove most of the QCD noise, having to deal finally with the semi-leptonic top 
pair decays which is the greatest challenge. As a matter of fact, di-jet pairs of light 
quarks from ti events naturally peak at M w ±, exactly as those from W ± H T do. As 
for the three jet systems, mispairings of b quarks with the wrong W ± have more severe 
effects on the signal than on the background, as one can intuitively expect from the 
production dynamics and as it can be appreciated in the middle plot of Fig. | (note 
the height of the peak for ti events, as compared to that of the W ± H T ones). 

Here we propose the following related and complementary cuts on invariant masses. 
For a start, we put ourselves in the favourable phenomenological position that the 
charged Higgs mass M#± is known, e.g., thanks to a previous detection of the light scalar 
Higgs h and to the measurement of its couplings. Under these circumstances, in order 
to enhance the W ± H Zf to ti rates, one can impose a cut on the invariant mass of the bb 
pair. Since in the signal both bottom jets originate from the H ± scalar and assuming 
that H ± —>■ bbW ± , the invariant mass squared M bb must be below J M^ ± — M^ ± (apart 
from finite width effects). The distribution from top pair events has a characteristic 
scale of 2mt and therefore, at low Mh±, it can be filtered out by setting a sufficiently 
low cut. The gain for the signal-to-background rate is large if H ± is reasonably (m t < 
M H ± < 2m t ) light (solid and dashed curves in Fig. ||, upper plot), whereas for heavy 
charged Higgs bosons, when M H ± is of the order 2m t or greater, the cut is not useful 
(dotted curve in Fig. |5], upper plot). However, one should note that this selection cut 
can be utilised successfully only when M#± is approximately known, and is of limited 
use even then, since it only removes the background from regions of phase space far 
away from reconstructing the charged Higgs mass peak. 

Similarly, one can impose a cut on the invariant mass of a bi pair, where b is the 
bottom jet which does not reproduce the top quark with the light di-jet pair (i.e., 
the one yielding the reconstructed m t further away from its actual value). Here, the 
selection works because for the ti background, if both top and are on the mass 
shell, one has M bf < ^mf — M^±. However, it should be noticed that the Higgs 

production mechanism can really push the Mm value beyond ymf — M^ ± only if M#± 
is large enough: see dashed and dotted curves in Fig. [| for M bi /j ~ 160 GeV. Failing 



this condition, the suppression against the signal itself can be quite large (e.g., an 
additional rejection factor of five for M H ± = 214 GeV at tan/? = 30). However, it 
turns out that such a constraint is definitely necessary to bring down the background 
rates to manageable levels (as it contributes with an additional factor of thirty or so to 
the suppression of top-antitop events), so that we employ it even at low M H ± values. 
In addition, we have made the following, more standard cuts: 

1. isolation of the two bottom jets from the light quark jets and from each other, 
as we tentatively set the azimuthal-pseudorapidity separation at AR bbb j > 0.7 
between themQ; 

2. for an isolated lepton, we impose the cut AR^/j > 0.4 between the lepton and 
all jets; 

3. the light quark jet pair mass Mjj within M^± ± 10 GeV; 

4. the light quark jet pair and a bottom quark jet combine at least once to a mass 
M bjj of m t ± 10 GeV; 

5. all one-particle pseudorapidities (of leptons, light and heavy quarks) are con- 
strained within a detector region of 2.5; 

6. the transverse momentum cut was set at 20 GeV for all jets and leptons, and for 
the missing transverse momentum as well. 

Fig. displays the total rates for signal and background after the above selection 
cuts have been enforced. Since the latter depends on the value of M H ± , the cross section 
for events of the type (|3J) is no longer constant when M#± ~ m t . However, given the 
weak dependence of Mh± on tan/5, the two top-antitop curves overlap in that mass 
range. Although the signal-to-background ratio has greatly improved, as compared 
to the initial situation in Fig. [1], for any combination of M#± and tan/3, this is still 
very small, at least one part in a hundred, so to presumably dash away any hopes of 
resolving the Higgs peak. 

In fact, to be realistic, one should expect a four jet mass resolution of no less than 
10 GeV, given the usual uncertainties in reconstructing parton directions and energies 
from multi-hadronic events. Therefore, in Fig. ^ we have binned the invariant masses 
of the bbjj system in signal and background using that value. If one does so, it is 
clear that the Higgs mass peaks at 214, 407 and 605 GeV (with Breit-Wigner width 
r#± of 1.2, 4.4 and 6.1 GeV, respectively), at tan/? = 30, are overwhelmed by the 
top-antitop events. Seemingly, even where the signal is more pronounced over the 
background (at large Mh±), the excess amounts to no more than 10% at the most in 
the central bin. This is probably too little, further considering, on the one hand, the 
aforementioned uncertainties (PDFs, if -factors, etc.) and, on the other hand, that 
the event rate is poor, only around 1 fb prior to the vertex tagging efficiency being 
applied. Needless to say, if one looks back at Fig. |7|, similar conclusions should be 
expected for all other combinations of Ma and tan/? considered here, finally recalling 
our systematic overestimate of the signal rates for low values of the latter. 

4 Note that we allow for the light quark jets to be arbitrarily close. 



4. Conclusions 



We believe that in the H T — > tb — > bbW T channel, heavy charged Higgs scalars of 
the Minimal Supersymmetric Standard Model produced in association with W ± gauge 
bosons cannot be resolved at the LHC, via semi-leptonic W + W~ decays, for Higgs 
masses in the range between 2m t and 600 GeV (those producible at observable rate), 
at neither low nor high tan j3, because of the presence of the irreducible background 
from top-antitop events. Furthermore, our results can safely be applied to a more 
general 2 Higgs Double Model too (where mass and coupling constraints in the Higgs 
sector can be relaxed), given the extremely poor significance of the W ± H :¥ rates over 
the ti ones. As for other hadronic collider environments, the prospects of detection at 
the Tevatron (Run II) are even more reduced, given the lower machine luminosity and 
since the production cross section of the signal is there about three orders of magnitude 
smaller than at the CERN accelerator (for detectable Higgs masses below 300 GeV or 
so), while the background only decreases by about two orders. 

We have reached these conclusions after performing a detailed signal-to-background 
analysis, based on matrix element calculations of elementary 2 — > 6 suprocesses, convo- 
luted with up-to-date parton distribution functions, and exploiting dedicated selection 
cuts, beyond the usual requirements in transverse momentum and pseudorapidity. Al- 
though we have confined ourselves to the the parton level only, wherein jets are identified 
with partons, we are however confident that hadronisation and detector effects will not 
modify our main results. 

Nonetheless, we would like to conclude this paper with a positive note. Charged 
Higgs production in association with W ±, s, via bb and gg fusion at hadron colliders, 
represents a novel mechanism, whose decay phenomenology is largely unknown and 
that should be investigated further, considering that the detection of this particle is 
not at all certain at the next generation of hadronic machines, especially in the heavy 
mass range. In this respect, we would like to advocate, for example, the consideration 
of non-Standard Model decay channels, involving squarks, sleptons and gauginos, which 
was beyond the intention of this study. 
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Figure 1: Total event rates for W^H 1 " (solid) and tt production (dashed) at the LHC, 
with no cuts implemented, using CTEQ(4L), as a function of M#± for tan/3 = 1.5 and 
tan p = 30. In the blow up figure, the W ± H T production rates for tan/5 = 1.5 including 
squark loop contributions, with rriq = 300 (dotted) and 500 (dot-dashed) GeV, divided 
by those obtained when the latter are neglected. For tan/5 = 30, squark contributions 
are negligible, so that the same ratios would visually coincide with one. 
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Figure 2: Differential spectra in lepton/light quark jet (top-left), b jet (top-right) and 
missing (bottom) transverse momentum. Here, tan (3 = 30 and the W ± H Zf rates are 
plotted for M H ± = 214 (solid), 407 (dashed) and 605 (dotted) GeV. As M H ± > m t +m b , 
ti events have no MSSM parameter dependence (dot-dashed). No cuts have been 
implemented. The PDFs CTEQ(4L) have been used. Distributions are normalised to 
unity. Note that, in the case of the lepton and light quark jet, the curves coincide. 
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Figure 3: Differential spectra in lepton/light quark jet (left) and 6 jet (right) pseudora- 
pidity. Here, tan (3 = 30 and the W ± H T rates are plotted for M H ± = 214 (solid), 407 
(dashed) and 605 (dotted) GeV. As M H ± > m t + mj, tt events have no MSSM param- 
eter dependence (dot-dashed). No cuts have been implemented. The PDFs CTEQ(4L) 
have been used. Distributions are normalised to unity. Note that, in the case of the 
lepton and light quark jet, the curves coincide. 
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Figure 4: Differential spectra in pseudorapidity-azimuth separation between the follow- 
ing pairs of particles in W^H^ and tt events: lepton/6 jet (top); lepton/light quark jet 
(bottom). Here, tan/3 = 30 and the W ± H Zf rates are plotted for M#± =214 (solid), 
407 (dashed) and 605 (dotted) GeV. As M H ± > m t + m b , tt events have no MSSM 
parameter dependence (dot-dashed). No cuts have been implemented. The PDFs 
CTEQ(4L) have been used. Distributions are normalised to unity. 
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Figure 5: Differential spectra in invariant mass for the following systems in W ± H T 
and tt events: bb (top); three jet (middle); four jet (bottom). Here, tan/3 = 30 and the 
W ± H T rates are plotted for M H ± = 214 (solid), 407 (dashed) and 605 (dotted) GeV. 
As M H ± > m t + mb, tt events have no MSSM parameter dependence (dot-dashed). No 
cuts have been implemented. The PDFs CTEQ(4L) have been used. Distributions are 
normalised to unity. 




Figure 6: Differential spectra in invariant mass for the b£/jet system in W ± H T and tt 
events. Here, tan/3 = 30 and the W ± H T rates are plotted for M H ± = 214 (solid), 407 
(dashed) and 605 (dotted) GeV. As M H ± > m t +mb, tt events have no MSSM parameter 
dependence (dot-dashed). No cuts have been implemented. The PDFs CTEQ(4L) have 
been used. Distributions are normalised to unity. 




Figure 7: Total event rates for W^H 1 " (solid) and ti production (dashed) at the LHC, 
after the selection cuts have been implemented, using CTEQ(4L), as a function of M#± 
for tan/? = 1.5 and tan/5 = 30. 
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Figure 8: Differential spectra in the four jet invariant mass in W^H^ and tt events 
for tan (3 = 30 and M H ± = 214 (solid), 407 (dashed) and 605 (dotted) GeV. Their sum 
(dot-dashed) is also reported in the three blow up frames below, compared against the 
ti rates from above. Even though M#± > m t + int,, ti events have a MSSM parameter 
dependence due to an M#± based constraint being implemented in the selection cuts. 
The PDFs CTEQ(4L) have been used. Distributions are normalised to total cross 
sections. 



